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Robust constraints on dark energy and gravity 
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ABSTRACT 

Galaxy clustering data provide a powerful probe of dark energy. We examine how the con- 
straints on the scaled expansion history of the universe, Xh(z) = H(z)s (with s denoting 
the sound horizon at the drag epoch), and the scaled angular diameter distance, Xd(z) = 
Da(z)/s, depend on the methods used to analyze the galaxy clustering data. We find that 
using the observed galaxy power spectrum, P° bs (k), Xh(z) and Xd(z) are measured more 
accurately and are significantly less correlated with each other, compared to using only the in- 
formation from the baryon acoustic oscillations (BAO) in P° bs (k). Using the {xh(z), Xd{z)} 

from P° bs (k) gives a DETF dark energy FoM approximately a factor of two larger than using 
the {xh(z), Xd(z)} from BAO only; this provides a robust conservative method to go beyond 
BAO only in extracting dark energy information from galaxy clustering data. 

We find that a Stage IV galaxy redshift survey, with 0.7 < z < 2 over 15,000 (deg) 2 , 

can measure {xh{z), Xd(z), f g {z)G{z)pQ 2 /s 4 } with high precision (where f g (z) and G(z) 
are the linear growth rate and factor of large scale structure respectively, and Pq is the dimen- 
sionless normalization of P° bs {k)), when redshift-space distortion information is included. 

The measurement of f g (z)G(z)pQ^ 2 /s i provides a powerful test of gravity, and significantly 
boosts the dark energy FoM when general relativity is assumed. 

Key words: cosmology: observations, distance scale, large-scale structure of universe 



1 INTRODUCTION 

The cosmic acceleration (i.e., dark e nergy) was discovered in 1998 
dRiess et al. 19981; IPerlmutter et alj|l9991) . and we are still in the 
dark about the nature of this mystery. We can hope to measure both 
the cosmic expansion history and the cosmic large scale structure 
growt h history accurately and precisely with galaxy clustering (see, 
e.g., IGuzzo et al.l ^20081): IWanj d2008ah) and weak le ns ing (see, 
e.g.. | Knox. Song. & Ty son (2006); Z hang etal] J2007h; iHeavens 
(2009)) data from a space mission such as Euclid jLaureiis et al. 
201 ljj, and differentiate the two possible explanations for the ob- 
served cosmic acceleration: a new energy component, or a modifi- 
cation of Einstein's theory of gravity^ 

Galaxy cl ustering has long been used as a cosmological 
probe (see, e.g., Hamilton (1998)). At present, the largest data set 
comes from the SPSS III Ba r yon Oscillatio n Spectroscop ic Survey 
(BOSS), see lAnderson et al.l 120121) and lReid et all J2012l) . 



* E-mail: wang@nhn.ou.edu 

1 http://www.euclid-ec.org 

2 Cluste rs of galaxies provide a c omplementary probe o f dark energ; 



Here we explore different analysis techniques for galaxy clus- 
tering data, in order to obtain robust constraints on dark energy and 
general relativity. It is important to study how the dark energy and 
gravity constraints depend on the assumptions we make about the 
information that can be extracted from galaxy redshift survey data. 

We present our methods in Sec. 2, results in Sec. 3, and sum- 
marize in Sec.4. 



2 METHODOLOGY 



2.1 The Fisher Matrix Formalism 



We use the Fisher matrix formali sm to study the parameter estima- 
tion using galaxy clustering data (Tegmark 1997; Seo & Eisenstein 
2003h. based on the approach developed in I Wand bOolT koOSa. 



2010al) : IWang et al] d2010h . In the limit where the length scale cor- 



see. e.g. jMaiumdar & Mohil l2004l) ; lManera & Motj l200el) ; lMotj (200 
ISartorisetalJfcOllh . 



responding to the survey volume is much larger than the scale of 
any features in the observed galaxy power spectrum P g (k), we 
can assume that the likelihood funct ion for the band powers of 
a gala xy redshift survey is Gaussian ( Feldman, Kaiser. & Peacock! 
1 19941) . with a measurement error in In P(k) that is proportional to 
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p4//(k)] 1 ' 2 , with the effective volume of the survey defined as 



dr' 



n(r)P g (k,n) 



n(r)P g (k,n) + l 



nP g (k,/i) 



nP g {k,n) + 1 



^survey •, 



(l) 



where the comoving number density n is assumed to only depend 
on the redshift (and constant in each redshift slice) for simplicity in 
the last part of the equation. 

In order to propagate the measurement error in In P g (k) into 
measurement er rors for the parameters pi, we use the Fisher matrix 
dTegmarkll 19971) 



F 



<91nP 9 (k) <91nP 9 (k) 



dpi 



dpj 



dk J 



2(2tt) 3 



(2) 



where pi are the parameters to be estimated from data, and the 
derivatives are evaluated at parameter values of the fiducial model. 
Note that the Fisher matrix Py is the inverse of the covariance ma- 
trix of the parameters pi if the pi are Gaussian distributed. 

The observed galaxy power spectrum can be reconstructed us- 
ing a particular reference cos mology, including the ef fects of bias 
and redshift-space distortions (Seo & Eisenstein 2003): 



P° hs {k r f ,k r { f ) 



'H(z) 



[D A (z)] 2 H(zy-f 
■\G{z)fP m {k) z=0 + P shot 



(3) 



where H(z) — a /a (with a denoting the cosmic scale factor) is 
the Hubble parameter, and Da(z) = r(z)/(l + z) is the angular 
diameter distance at z, with the comoving distance r(z) given by 



r(z) = c\Q k \ 1/2 si: 



I fij 



,1/2 



dz 



H(z') 



(4) 



where sinn(s) = sin(:r), x, sinh(:r) for fife < 0, fife = 0, and 
fife > respectively. The bias between galaxy and matter distribu- 
tions is denoted by b(z). Th e linear reds hift-space distortion param- 
eter /3(z) — f g (z)/b(z) dKaisenl 19871) . where f g (z) is the linear 
growth rate; it is related to the linear growth factor G(z) (normal- 
ized such that G(0) = 1) as follows 



/»(*) 



d In G(z) 



(5) 



d In a 

Note that fi = k ■ r/fc, with f denoting the unit vector along 
the line of sight; k is the wavevector with jkj = k. Hence /i 2 = 

k\/k 2 = fcjj/(fci + k\), where 



fc|| 



k ■ f = kfi 



yjk* - *J = ky/T- 



(6) 



The values in the reference cosmology are denoted by the subscript 
"ref", while those in the true cosmology have no subscript. Note 
that 

k r l s = k ± D A (z)/D A (z) ref , k\ ef = k^H{ Z ) rBf /H{z). (7) 

Eq.lO characterizes the dependence of the observed galaxy 
power spectrum on H(z) and Da(z) due to BAO, as well as the 
sensitivity of a galaxy r edshift survey to the linear redshift-space 
distortion parameter f3 jKaiseri [198% . The linear matter power 
spectrum at z — is given by 



P m (fc) 2 =o = Pofc" s T 2 (fc), 



(8) 



where T(k) denotes the matter transfer function. 

The measurement of f g (z) given /3(z) = f g (z) /b(z) requires 
an additional measurement of the b ias b(z), which could be ob- 
tained from the galaxy bispectrum dMatarrese. Verde. & Heavens! 
1 19971 ; I Verde et al]|2002l) . Alternatively, we can rewrite Eq.([3j as 



PS bB {k r l ! ,k r ^ ! ) 

p° bs (k™ f ,k: ef )/(h- 1 Mpc) 3 



\D A (zy-f} 2 H{z) 



[D A (z)] 2 H(z)™f \Mpc~ 1 
■ [o g (z) + f g (z)a m (z) /i 2 ] 2 + P shot 

where we have defined 

a g (z) = b(z) G(z) P 1/2 
a m (z) = G(z)P 1/2 , 



T 2 (fc) 



(9) 



(10) 



The dimensionless power spectrum normalization constant Po is 
just Po in Eq.l[8]l in appropriate units: 



Po 



Po 



- 2 



(11) 



(Mpc/ft) s (Mpc)™ s I h n °' 

The second part of Eg. ill It is relevant if as is used to normalize the 
power spectrum. Note that 



Jo 



dk — T 2 (k- hMpc' 1 ) 

2tt z 



3ji(8fc) 



(12) 



where k = k/[hMpc 1 ], and ji(kr) is spherical Bessel function. 
Note that Io = Io(u} m , Wi>, n s , h). Since ku and k± scale as H(z) 
and 1/Da(z) respectively, P g bs (k) in Eq.([9j does not depend on 
h. 

Eq.([9]( is analogous to the approach of ISong & Percivall 

(2009), who proposed the use of f g (z)as(z) to probe growth 
of large scale structure. The difference is that Eq.((9} uses 
f g (z)am(z) = f g (z)G(z)Pg , which does not introduce an ex- 
plicit dependence on h (as in the case of using f g (z)(r&(z)). 

The uncertainty in redshift measurements is included by mul- 
tiplying P g (k) with the damping factor, e~ k M " r , due to redshift 
uncertainties, with 



dr 

o r = — a z 
az 



(13) 



where r is the comoving distance. Note that the damping factor 
should be held constant when taking derivatives of P g (k). 

2.2 P(k) Method 

Including the nonlin ear effects explicitly, we can write 
dSeo & Eisensteinll2007h 



dP g {k,n\ z ) _ dP l g in {k,n\z) 



exp 



(14) 



dpi dpi 

The damping is applied to derivatives of P g (k), rather than P g (k), 
to ensure that no information is extracted from the damping itself. 
Eq.l(2j becomes 



Fij — Vsurvey 



1 rkma * 91nP^ n (fc )Al ) d\nP l g m (k,y) 



1 Jk„ 



dp 3 



nP l g ln (k,v) 



nP l g in (k,n) + 1 



2 S 2 ; 27Tfc 2 dfc 

2(2tt) 3 ' 



(15) 
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The linear galaxy power spectrum P g 
The nonlinear damping scale 



'(fc, y\z) is given by Eq.l[9]l. 



d\nx/dp2 — fl 



(24) 



En 



(1 - fj, )E ± + ;u 
E±(l + / 9 ) 



Ei = 12 Ah 



F V0.9 



8.355 fc _1 Mpc 



0.758 G(z)p NL 
■ G{z)p NL . 



(16) 



The parameter pnl indicates the remaining level of nonlinear- 
ity in the data; with pnl = 0.5 (50% nonlinearity) as the 
best case, and Pnl = 1 (100% nonlinearity) as the worst case 
dSeo & Eisensteinll2007l). F or a fiducial model based on WMAP3 
results dSpergel et : alj|2007h (fl m = 0.24, h = 0.73, Ov = 0.76, 
n k = 0, Q. b h 2 = 0.0223, r = 0.09, n s = 0.95, T/ S = 0), 
A = 0.5817, P .2 = 2710ct| i9 dSeo & Eisensteinll2007h . 

In the P(k) method, the full set of parameters that de- 
scribe the observed P g (k) in each redshift slice are: {In H(zi), 
InDA(zi), h\[f g (zi)a m (zi)], \n<7 g {zi), Pl hot \ Mm, w&, n s }, 
where Lu m = Q m h 2 , and u) b = Qth 2 . We marginalize 
over {In a g ( Zi), Pg hot } in each redshift slice, to obtain a Fisher 
matrix for {lnH(zi), In Di(zj), ln[f g (zi)(T m (zi)]; w m , Ub, n 3 }. 
This full Fisher matrix, or a smaller set marginalized over various 
parameters, is projected into the standard set of cosmological pa- 
rameters {wo,w a , &x , , u m , w& , « s , In As } ■ There are four dif- 
ferent ways of utilizing the information from P(k) (see Sec l3.3t . 

2.3 BAO Only Method 

The power of galaxy clustering as a dark energy probe was first 
recognized vi a studies of baryon acoustic oscillations (BAO) as a 
standard ruler dBlake & GlazebrookfeOOllSeo & Eisenstei j2003h . 
The BAO only method essentially approximates dP l g n (k, fi)/dpi 
with dPl in (k,^\z)/dpi in the derivatives in Eq.fTSJ, with the 
po wer spectrum that contai ns baryonic features, Pl' n (k, fi), given 
by dSeo & EisensteiiJl2007l) 



The BAO only method gives {xh(z), x d (z)} that are correlated at 
the level of ~41% with each other, but are uncorrelated for different 
redshift slices by construction. 



2.4 Assumptions and Priors 

We u se the fiducial model adopted by the FoMSWG dAlbrecht et al.l 
l2009h . with ui m = Q, m h 2 = 0.1326, uj b = Q. h h 2 = 0.0227, 
h = 0.719, Q k = 0, w = -1.0, n s = 0.963, and a 8 = 0.798. No 
priors are used in deriving {xh(z),Xd{z), f g (z)a m (z)/s 4 }, which 
provide model-independent constraints on the cosmic expansion 
history and the growth rate of cosmic large scale structure. These 
allow the detection of dark energy evolution, and the differentia- 
tion between an unknown energy component and modified gravity 
as the causes for the observed cosmic acceleration. 

In order t o derive dark energy figure of merit (FoM), as defined 
by the DETF dAlbrecht et al.ll200q) , we project our Fisher matrices 
into the standard set of dark energy and cosmological parameters: 
{w Q ,w a ,QxMk,u m ,oj b ,n a AnA s }. To include Planck priors0 
we convert the Planck Fisher matrix for 44 parameters (including 
36 parameters that parametrize the dark energy equation of state in 
redshift bins) from the FoMSWG into a Planck Fisher matrix for 
this set of dark energy and cosmological parameters. 

We present all our results for StagelV+BOSS spectroscopic 
galaxy redshift surveys. The Stage IV galaxy redshift survey is as- 
sumed to cover 15,000 (deg) 2 , with Ha flux limit of 3 x 10~ 16 
ergs _1 cm -2 , an efficiency of e = 0.50, a redshift range of 
0.7 < z < 2.05, and a redshift accur acy of a z /(l + z) = 0.001. 
The galaxy number density is given by Geach et al. and the 

galaxy bias function is given bv lOrsi et alj d20101) . This is similar 
to the baseline of the Euclid galaxy redshift survey dLaureiis et al.l 
l201ll) . The BOSS survey is assumed to cover 10,000 (deg) 2 , a red- 
shift range of 0.1 < z < 0.7, with a fixed galaxy number density 
of n = 3 x 10~ 4 /i 3 Mpc~ 3 , and a fixed linear bias of b = 1.7. 



Pb"(k,n\z) 



^A P^(k .2,tJ,\z) S -^-exp [-(fcE,) 1 ' 4 ] , (17) 



where P l g in {k, n\z) is the linear galaxy power spectrum, and the 
Silk damping scale E s = 8.38 h~ Mpc. We have defined 



/, 2 2 . , 2 2\l/2 

[k ± s ± + fen sii ) 



fco.2 = 0.2 h Mpc" 1 

x = 
Defining 

Pi = lns^ 1 = \n(D A /s) = \n(x h ) 



(18) 
(19) 



3 RESULTS 

3.1 Measurement of H(z) and Da(z) 

The H(z) and Da(z) observables that correspond to the BAO scale 



(20) 

p 2 = lns|| =bx(sH) = ln(x d ), (21) 

substituting Eq.dl7b into Eq.dl5b. and making the approximation of 
cos 2 x ~ 1/2, we find 

1 r k„ 

A 2 I 

zj — v survey-' 



Xh(z) 
Xd(z) 



H(z)s 
D A (z)/s 



(25) 



Fa ~ VsurveyAQ I d(j,fi(n)fj(n) / dfefc 2 



(1 



P l ™{k\z = 0) 



+ 



Pik n {k . 2 \z = Q) ' nP^ n (ko.2,fJ-\z)e~ k2 ^"r 
•exp [-2(feE s ) 1 ' 4 -fe 2 E 2 ! ] , 

where P""(fe . 2 , jj,\z) is given by Eq.l|9j with fc = fco.2. 



The functions /; (p) are given by 



dlnx/dpi 



(22) 



(23) 



where s is the sou nd horizon scale at the drag epoch 
dHu&Surivamalll996l) . 

Fig[T] shows the measurement precision of Xh{z) and Xd(z) 
for StagelV+BOSS. The top panel shows the percentage errors on 
Xh(z) and Xd(z), the bottom panel shows the normalized corre- 
lation coefficient between them. The thick solid and dashed lines 
represent the measurement precision of Xh(z) and x d (z) from the 
P(k) method, marginalized over all other parameters. The thin dot- 
ted and dot-dashed lines represent the measurement of Xh{z) and 
Xd(z) from the BAO only method. 

Note that the Xh{z) and Xd(z) measured using P(k) are 



3 For a genera l and robust method for including Planck priors, see 
Mukhe riee etaf! 120081) . 
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— i — i — | — i — , — i — , — 

■ Alnx h (z), P(k) method 
Alnx„(z), BAO method 

■ Alnx d (z). P(k) method 
Alnx d (z), BAO method 




0.2h/Mpc, p Nl =50% 



StagelV + BOSS 



-0.5 



P(k) 



=0.2h/Mpc, p„ L =50% 



Figure 1. Precision of Xh (z) = H(z)s and Xd(z) = Da(z)/s expected 
from StagelV+BOSS. The top panel shows the percentage errors on x n (z) 
and x d (z) per Az = 0.1 redshift bin, the bottom panel shows the normal- 
ized correlation coefficient between x/ l (z) and Xd(z). 



only weakly correlated. Since x n (z) and Xd(z) represent indepen- 
dent degrees of freedom in a galaxy redshift survey, they should 
not be strongly correla ted. This is consistent with the findings of 
IChuang & wingl fcOllh from their analysis of SDSS LRG data. 

The Xh(z) and Xd(z) from using the BAO method are corre- 
lated with a normalized correlation coefficient of r ~ 0.41. They 
are positively correlated by construction: In both P(k) and BAO 
only methods, the Fisher matrix element for {xh(z), Xd(z)} from 
the same redshift slice is negative. In the P(k) method, the 11,(2) 
and Xd(z) from different redshift slices are correlated through the 
cosmological parameters {o; m , lu b , n s } that are measured using in- 
formation from all the redshift slices. When the cosmological pa- 
rameters are marginalized over, the dependence on these parame- 
ters remain as a weak correlation between Xh(z) and Xd(z). In the 
BAO method, the Xh(z) and x d (z) from different redshift slices 
are uncorrelated by construction. Inverting the 2x2 Fisher matrix 
of x h(z) and Xd(z) in each redshift slice leads to positive and sig- 
nificant correlation between xu(z) and Xd{z). 

Finally, note that in using the P(k) method to forecast dark 
energy constraints, two different methods have been used to ac- 
count for nonlinear effects: 

(1) Setting kmax = ir/(2R), with R given by requiring that <J 2 (R) 
is small (e.g., & 2 (R) = 0.25), and imposing a uniform upper limit 
cutoff, e.g., kmax ^ 0.2ft/Mpc. Note that in this casepjvL = in 
Eq.dl6|>; the nonlinear effects are minimized by imposing a mini- 
mum length scale that increases at lower redshift. 

(2) Setting k max to a fixed value, and account for nonlinear effects 
through the exponential damping term in Eq.dl6l>, e.g., pnl = 0.5. 



With a suitable choice of cr 2 (7?) in (1) and k max in (2), 
these two methods of accounting for nonlinear effects give the 
same DETF dark energy FoM. For StagelV+BOSS, for the P(k) 
only method (no priors and marginalizing over growth informa- 
tion), (1) with a 2 (R) = 0.28 and k max ^ 0.2/i/Mpc gives 
FoM=49.9, while (2) with p NL = 0.5 and k max = 0.2/i/Mpc 
gives FoM=49.6. These two cases give very similar uncertainties 
on Xh(z) and x d (z). 

Since these two nonlinear cutoff methods are very similar, we 
have chosen to use cutoff method (2) in the rest of this paper, since 
it is smooth with k, and is the approach used in the BAO only 
method. 



3.2 Growth Rate Measurements 



ISong & Percivall J2009I) showed that assuming a linear bias be- 
tween galaxy and matter distributions, we can use f g (z)as(z) to 
probe gravity without additional assumptions. We use a similar ap- 
proach, but use f g (z)a m (z) = f g (z)G(z)py 2 to avoid intro- 
ducing an explicit dependence on h through as- We find that the 
fg(z)tJm(z) measurements from the P(k) method are highly cor- 
related with the u m measurement, which makes the uncertainties 
on f g (z)a m (z) much larger than that of /3(z). Fortunately, we are 
able to find a scaled measurement of f g (z)a m (z), 

f g (z)a m (z) _ f g (z)G(z)P^ 2 



f g {z)a m (z) = 



(26) 



that is nearly uncorrelated with uj m , and has an uncertainty that ap- 
proaches that of /3(z), see top panel of Fig[2] The precision of both 
f g (z)a m {z)/ s 4 and /3(z) are insensitive to the choice of k max - 

To make sense of the scaling in Eq.J26l>. note that the observed 
power spectrum depends on uj m only through s and T(k), as fol- 
lows (see Eq.(9)) for P shot = 0: 



Xh(z) 



1 r 2 

— ■ [0-9(2) + f g (z)a m (z) /J, 



(ks 



- lO g {z) + f g (z)a m (z) v 2 ] 2 T 2 (k)s-^ +3 \ 
Note that at the peak of P m (k\ z = 0) = P k ns T 2 (k), k = k p , 
91nT(fc) 



-A + 2 



Ok 



= 0, 



(28) 



where k = fc/Mpc . T(k) depends only on 

k 



(-dm 



hMpc- 1 

where 62.7 = T C mb/2.TK, T 
Thus at k = k p , 

dlnT(fc) ^ 1 dlnT 

dk U m dq 

and we find 

dlnT dlnT dq 



(29) 

Qmh jEisenstein &~Hu]|l998h . 

(30) 



2oj m 



dlnT 
dq 



(3D 



where we have used Eqs.ll28t and OOt . Using the approximate for- 
mula for s from Eisenstein & Hu (1998), 
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..,44.5]n(9.83/^)„ r . 



3/4 



1 + 10W, 



we find that at the peak of P(k), 



dlnT 



du) n 



— ln(9.83/w m ) 

£ OUJm 

„ 91ns 



<9cj„ 



(32) 

(33) 
(34) 



We have assumed ui m and n s close to our fiducial values of lj„ 
0.1326 and n 3 = 0.963 in obtaining Eq.CT. 
If we define the scaled parameters 



a 9 («) _ 6(^)G( 2 )Po 



1/2 



(35) 



fg(z)<T m (z) = 

we find 



f g (z)Om(z) _ f 9 {z)G{z)P' 



1/2 



(2) + f g (z)a m (z)n -{ks) r 



s T 2 (fc).(36) 



In the new set of parameters, {xh(zi), Xd(zi), fg{zi)cr m {zi), 
a g (zi), Pl hot \ Lu m , cot, n s }, the dependence of Pg bs on ui m only 
comes through the combination of s 5 ~" 3 T 2 (k) ~ [s 2 r(fc)] 2 , 
which is only very weakly dependent on oj m (see Eg. 1341 ). Thus the 
dependence of Pg bs on io m is effectively removed or absorbed via 
the scaling of parameters in Eq.l|35l>, leading to measurements on 
fg(z)(Tm(z) that are essentially uncorrelated with u) m , and greatly 
improved in precision over that of f g (z)a m (z). This is as expected, 
since the measurements of f g (z)a m (z) are strongly correlated with 
that of Lu m (i.e., P(k) shape). 

The bottom panel of Fig|2] shows the uncertainties on the 
growth rate powerlaw index 7 for StagelV+BOSS, with and 
without Planck priors. Note that 7 is defined by parametriz- 
ing the growth rate as a powerlaw dWang & S teinhardt 1 19981 ; 
iLue. Scoccimarro . & Starkman 2004), 



/,(*) = [fi m (a)P 



(37) 



where fi OT (a) = 8ivGpm{a)/(3H 2 ). The solid lines in the 
bottom panel of Fig|2] show the precision on 7 using only 
the {xh(z), Xd{z), f g (z)a m (z)/s 4 } measured from P(k) and 
marginalized over all other parameters. The dashed lines show the 
precision on 7 when the full P(k) is used, including the growth 
information (i.e., the "P(k) + f g " method). 



3.3 Dark Energy Figure of Merit 

To calculate the DETF dark energy FoM jAlbrechtetaT] [2006), 
FoM= l/- v /det[Cov(ui , w a )] (Wang 20083), we need to project 
our large set of measured parameters into the standard set of 

{»(), W a , fix, fife, CJm,^, Us, In As}. 

The BAO only method gives measurement of 
{xh(zi),Xd(zi)} from each redshift slice. The Fisher matrix 
for these measurements are then projected into the Fisher ma- 
trix for {wo, w a , fix, fife, w OT , cjf,}. Because the dependence 
on {u)m,u } only comes through s, ui m and w& are perfectly 
degenerate if no priors are added. This can be shown explicitly by 
computing the submatrix for {w m , Mb} in the Fisher matrix, which 
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CD 
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Figure 2. Top: uncertainties on f g (z)G(z)P 1 ^ 2 / s 4 and j3(z) for 
StagelV+BOSS per Az = 0.1 redshift bin. Bottom: uncertainties on the 
growth rate powerlaw index 7 for StagelV+BOSS, with and without Planck 
priors. 



is proportional to 



'Sin 



' d In 5 



')( 



8 In 



n x h \ I d In 



(38) 



the determinant of this submatrix is zero, thus the determinant of 
the entire Fisher matrix for {wo, w a , fix, fife, Wm, oJb} is zero. It 
can be shown that the combination determined by the BAO only 
method is 

-1 



/ d In x h \ f dlnxh \ 
V dcob I V dui m J 



dbJb 

/d\ns\/ 1 + 91ns\ _1 
V du>b J \ 2uj m dujm. J 
It can be shown explicitly that the Fisher matrix for 
{wo, w a , fix, fife, OJm} is exactly the same as the 
{wo, w a , fix, fife, oj m } submatrix of the original Fisher ma- 
trix for {wq, w a , fix, Sli,w m ,wi). Thus to compute the FoM for 
BAO only, one only needs to drop the Fisher matrix elements for 
ujb, then invert the resultant Fisher matrix to obtain the covariance 
matrix. Note that the Fisher matrix for {wo,w a , fix, fife, OJm, w&} 
should be used when combining with Planck priors. 

There are four different ways that we can extract dark energy 
information from the P(k) method (in the order of increasing 
information content): 

(1) {xh(z), Xd(z)} from P(k): Project the Fisher matrix 
for {\nH(zi), \nD A {zi), ln[/ 9 (z l )cr m (2 i )]; u m , uj b , n s } 
into that of {Inxh(zi), \nXd{zi), ln[/g(2i)cr m (zi)/s 4 ]; 
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cj m , cjf,, n s } (see Sec l3.2b . then marginalize over 
{ln[fg(z i )a rn (z i )/s 4 ];Lu m ,Lo b ,n s }, and project the Fisher 
matrix for {In Xh(zi), IriXd(zi)} into the standard set of cosmo- 
logical parameters. 

(2) {xh(z),Xd(z), fg(z)a m (z)/s 4 } from P(k): Project the Fisher 
matrix for {InH(zi), ]nDA{Zi), ln[f g (zi)a m (zi)]; u) m , Wb, n a } 
into that of {lnxh(zi), \nx d (zi), ln[f g (zi)a m (zi) / s 4 ]; cj m , ui b , 
n 3 } (see Sec |3.2t , then marginalize over {cj m ,u>t,, n 3 }, and project 
the Fisher matrix for {lnx h (zi),lnx d (zi),ln[fg(zi)a m (zi)/s 4 ]} 
into the standard set of cosmological parameters. 

(3) P(k) marginalized over f g : Marginalize over ]n[f g (zi)a m (zi)] 
to obtain the Fisher matrix for {In H(zi), \nDA(zi)\u m ,u)b, n s }, 
and project it into the standard set of cosmological parameters. 

(4) P{k) + f g : Project the Fisher matrix for {In H(zi), In D A (zi), 
ln[fg(zi)a m (zi)]; u m , U)i, n s } into the standard set of cosmologi- 
cal parameters. 

Fig[3]shows the DETF dark energy FoM for StagelV+BOSS, 
without (top panel) and with (bottom panel) Planck priors, as a 
function of the nonlinear cutoff k max . The four methods of us- 
ing P(k) described above, as well as the BAO only method, 
are shown. Note that the FoMs from the three most conser- 
vative methods, BAO only, {xh(z), x d(z)} from P(k), and 
{xh(z),Xd(z), f g (z)am(z)/s 4 } from P(k), are insensitive to the 
increase of k max for k max > 0.3/i/Mpc. We have not included the 
nonlinearity in the RSD due to peculiar velocities here for simplic- 
ity. Adding a peculiar velocity of 300 km/s is equivalent to adding 
0.001 in quadrature to the redshift dispersion a z — 0.001(1 + z)\ 
this has a negligible effect on the FoM for StagelV+BOSS, since 
the FoM is most sensitive to assumptions about the Stage IV survey, 
which is at z 0.7. 

The most conservative of the P(k) approaches, using 
{xh(z),Xd(z)} measured from P(k) and marginalized over all 
other parameters, gives a dark energy FoM about a factor of two 
larger than that of the BAO only method, with or without Planck 
priors. This provides a robust conservative method to go beyond 
BAO only in extracting dark energy information from galaxy clus- 
tering data. 

It is interesting to note that {xh(z), Xd(z)} from P(k) (solid 
line) gives similar dark energy FoM to that of the full P(k) 
marginalized over growth information (dotted line), when Planck 
priors are included. Similarly, {xf l (z),Xd(z), f g (z)a m (z)/s 4 } 
gives dark energy FoM close to that of the full P(k) with growth 
information included, when Planck priors are added. 



3.4 Comparison With Previous Work 

This work has the most overlap with lWang et al.l bold) , which ex- 
plored the optimization of a spa ce-based galaxy re dshift survey. 
The differences of this work from lWang et al.l ( feoioh are: 

(1) This work presents a new conservative approach to extract dark 
energy constraints from galaxy clustering data: the use of only 
the H(z)s and Da(z)/s measurements from the observed galaxy 
power spectrum, Pg bs (k), to probe dark energy. This bridges the 
methods using P° bs {k) and the BAO method (which uses H(z)s 
and Da (z) / s measurements from fitting the BAO peaks). 

(2) This work presents a new combination of growth information, 
f g (z)G(z)/s 4 , that can be measured nearly as precisely as the 
linear redshift-space distortion parameter /3 (see Fig.2), but can 
be used to probe the growth history of cosmic la r ge sca le struc- 
ture without assuming a bias model. Wa ng et ail d2010h did not 
study growth constraints explicitly; they either marginalized over 
the growth rate information, or assumed that gravity is described 
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Figure 3. Dark energy FoM for StagelV+BOSS, without (top panel) and 
with (bottom panel) Planck priors, as a function of the nonlinear cutoff 

h 

h in a x ■ 



by general relativity. 

(3) This work focuses on model-independent constraints of dark en- 
ergy and gravity in terms of H(z)s, Da(z)/s, and f Q ( z)G(z) / s 4 
measured in Az = 0.1 redshift bins. IWangetal .1 (2010) focuses on 
the conventional dark energy model with d ark energy equation of 
state given by wx(z) = wo + w a (l — a) dChevallier & Polarskil 
l200ll) . and dark energy density function X(z) = px(z) / px{0) 
parametrized by its value at z — 2/3, 4/3, and 2. 

The methodology developed in this work differs from what 
is currently used in analyzing galaxy clustering data. This work 
proposes the simultaneous measurement of H(z)s, Da(z)/s, and 
f g (z)G(z)/s 4 from galaxy clustering data without imposing any 
priors. Because of the limited volume probed by current data, no 
simultaneous measurements of H(z), Da(z), and f g (z) have been 
made without imposing strong priors on cosmological parame- 
ters. The first si multaneous measureme nts of H[z)s and Da{z)/s 
were made by IChuang & Wand d201ll) at z e f / = 0.35 using 
SPSS DR7 LRG d ata; they marginalized over growth information. 
iBlake et al] J201 lb measured f g {z)a%{z) at seve ral redshift s while 
fixing the background cosmology. Most recently. iReid et ai]d2012h 
published the first simultaneous measurement of Da(z)H(z), 
D a (z) 2 /H(z), and f g (z)a$(z) at z eff = 0.57 using BOSS data, 
assuming WMAP7 priors. 

This work presents forecasts of the precision of the most gen- 
eral measurements of cosmic expansion history (via H(z) and 
Da{z)) and gravity (via f g (z)G(z)) that can be made from a 
Stage IV galaxy redshift survey. These will be the most model- 
independent results on probing dark energy and probing gravity 
from such a survey. 
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4 SUMMARY AND DISCUSSION 
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We have examined how the constraints on the scaled expansion 
history of the universe, n(z) = H(z)s, and the scaled angu- 
lar diameter distance, Xd(z) = Da(z)/s, depend on the meth- 
ods used to analyze the galaxy clustering data. We find that using 
the observed galaxy power spectrum, Pg bs (k), xh{z) and Xd{z) 
are measured more accurately and are significantly less correlated 
with each other, compared to using only the information from the 
baryon acoustic oscillations (BAO) in Pg bs (k) (see Fig.l). Using 
the {x h (z), x d (z)} from P° bs (k) gives a DETF dark energy FoM 
approximately a factor of two larger than using the {xh(z), Xd(z)} 
from BAO only (see Fig. 3); this provides a robust conservative 
method to go beyond BAO only in extracting dark energy in- 
formation from galaxy clustering data. This is encouraging since 
IChuang & Wang! J201 lh found that {x h (z), x d (z)} from SDSS 
galaxy clustering data are not sensitive to systematic uncertainties. 

Furthermore, we find that if the redshift-space distortion 
information contained in P° bs (k) is used, we can measure 

!/2 ( 



{x h {z),Xd(z), fg(z)G(zji o 



Pn'^/s } with high precision from a 



Stage IV galaxy redshift survey with 0.7 

n2 



< z < 2 over 15,000 
(deg)^ (see Figs.l and 2), where f g (z) and G(z) are linear 
growth rate and growth factor of large scale structure respectively, 
and Pa denotes the dimensionless normalization of P^ n (k\z = 
0) . Adding f g (z)G(z)pQ /2 / s 4 to {x h (z), x d (z)} significantly 
boosts the dark energy FoM, compared to using {xf l (z),Xd(z)} 
only, or using P° bs {k) marginalized over the growth informa- 
tion, assuming that gravity is not modified (see Fig.3). Alter- 
natively, f g (z)G{z)P^ 2 / s A provides a powerful test of grav- 
ity, as dark energy and modified gravi ty models tha t give iden- 
tical H(z) likely give different f g (z) (Wang 2008a). Measuring 
{xh(z), Xd(z), f g (z)G(z)P^ 2 / s 4 } simultaneously allows us to 
probe gravity without fixing the background cosmological model. 
We will be adopting this approach to analyze simulated and real 
galaxy redshift catalogs in future work. 

We have developed a conservative approach to analyzing 
galaxy clustering data that should be insensitive to systematic un- 
certainties, if only data on quasi-linear scales are used (k ma x < 
0.2/i/Mpc). Since the dark energy FoM (see Fig.3) and the gravity 
constraints (see lower panel of Fig.2) are insensitive to the inclu- 
sion of smaller scale information at k max > 0.3/i/Mpc, our results 
are likely robust indicators of how well a Stage IV galaxy redshift 
survey can probe dark energy and constrain gravity. 

In analyzing real data, the systematic effects (bias between 
luminous matter and matter distributions, nonlinear effects, and 
redshift-space distortions^ will ultimately need t o be modeled 
in det a il and redu c ed w h ere possible (see, e . g., Percival et al.l 
d2010l) : lBlake et all d201ll) : |Padmanabhan et alj J2012h ). This will 
require cosmological N-body simulations that include galax- 
ies, either by inco rporat i ng ph y sical models of ga laxy for- 
mation (see, e.g., Baugh (2006); Angul o et all j2008h ). or us- 
ing halo occupation distributions (HOD ) measur e d fro m the 
largest available data sets (see, e.g., IZheng et all J2009h . and 
http://lss.phy.vanderbilt.edu/lasdamas/overview.htmll. We can ex- 
pect a Stage IV galaxy redshift survey to play a critical role in ad- 
vancing our understanding of cosmic acceleration within the next 
decade. 



4 See, e. g iBlake & Glazebrookl j2003l);ISeo & Eisensteij <2003l) . For re- 
el^ 



views, see 
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